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RESUMEN 

Se realiza el análisis sísmico de un edifico, ubicado en la ciudad de Manta 
con un modelo que considera base empotrada, aplicando el método espectral se 
obtienen desplazamientos y fuerzas en el centro de masas de la estructura, para lo 
cual se emplea el espectro de Manta en componentes E-W y N-S. Luego se 
presenta la teoría para encontrar la respuesta dinámica espectral de cualquier 
pórtico que conforma la estructura. Por otra parte se utiliza la metodología de la 
ASCE-10 para obtener un modelo simplificado de un grado de libertad equivalente 
y la guía NEHRP-12 para considerar interacción suelo estructura, además, se 
detalla el cálculo de los factores que se tomaron en cuenta para el modelamiento 
de una estructura con vigas de cimentación; y, se comparan los resultados de los 
distintos modelos propuestos. Por último, en base a la deriva de piso se explica el 
daño registrado durante el terremoto del 16 de abril de 201 6, de magnitud 7.8. 

Palabras Clave: Método espectral, modelo simplificado de un grado de libertad 
equivalente, interacción suelo estructura. 


ABSTRACT 

Seismic analysis of a building located on Manta city with a model that 
considers embedded base is made, applying the spectral method displacements 
and torces in the center of masses of the structure are obtained, for which the 
spectrum of Manta is used in components E-W and N-S. Then the theory is 
presented to find the spectral dynamic response of any frame that conforms the 
structure. On the other hand the methodology of ASCE-10 is used to obtain a 
simplified model of an equivalent degree of freedom and the NEHRP-12 guide to 
considerer soil structure ¡nteraction, in addition, is detail the calcularon of factors 
that were taken into account for the modeling of a structure with strip footings; and, 
the results of the different proposed models are compared. Finally, based on the 
story drift the damage recorded during the earthquake of April 16, 2016 of 
magnitude 7.8 is explained. 
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KeyWords: Spectral method, simplified model of an equivalent degree of freedom, 
soil structure interaction. 


1. INTRODUCCIÓN 

Durante el siglo XX la costa norte de Ecuador soportó una serie de sismos 
de alta magnitud, el gran terremoto en 1906 (Mw 8.8) ocasionó la ruptura de 
aproximadamente 500 km de longitud de la placa principal; en las siguientes 
décadas en la misma zona se rompieron diversos sub-segmentos a causa de 
sismos suscitados de manera secuencial de sur a norte en: 1942 (Mw 7.8), 1958 
(Mw 7.7) y 1979 (Mw8.2). 

El terremoto de Pedernales en 2016 (Mw 7.8) supone una nueva etapa de 
eventos sísmicos, debido a que se encuentra en el mismo lugar del segmento ya 
roto en el terremoto de 1942 (Mw 7.8), estos acontecimientos fueron el resultado 
del empuje de fallas poco profundas en el límite de las placas de Nazca y 
Sudamérica; en el lugar del terremoto, la placa de Nazca subduce hacia el este 
por debajo de la placa Sudamericana a una velocidad de 47 mm / año. 

La figura 1 se refiere a los incidentes sísmicos mencionados. La línea 
negra con triángulos señala la fosa de subducción; las estrellas negras muestran 
los lugares de los epicentros de los grandes terremotos del siglo XX; los círculos 
blancos gruesos delimitan el área de liberación de energía de los terremotos; la 
línea blanca paralela a la costa ecuatoriana muestra la extensión de ruptura 
generada por el gran terremoto de 1906 (Mw 8.8); la estrella amarilla marca el 
epicentro del terremoto del 16 de abril de 2016 en Pedernales. 



Longitude 


Figura 1 Zona de eventos sísmicos de Ecuador. 
Fuente: Hayes et al. (2016) 


El sismo de Pedernales rompió un segmento de cerca de 100 km de 
longitud por 40 km de ancho de la placa principal en profundidades entre 15 y 30 
km, esto dio lugar a un deslizamiento sísmico de norte-sur de uno a seis metros lo 
que ocasionó como resultado la ruptura de dos asperezas adyacentes distintas. La 
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segunda aspereza de 40x40 km2 se rompió en un lapso de tiempo de 10 s con un 
deslizamiento superior a 5 m y generó más de la mitad de la liberación total de la 
energía sísmica, esta parte de la ruptura situada a 20 km de profundidad de la 
zona poblada de la costa oeste hace relación al extenso daño sufrido de 
Pedernales a Canoa. (Nocquet, 2016). 


En la figura 2 se presenta la evolución de la ruptura en un espacio de 
tiempo, a) acumulación de la energia cada 6 s, las líneas rojas son los contornos 
de profundidad de la subducción cada 20 km; b) función tiempo versus energía 
liberada; c) tres fases principales de ruptura en los intervalos de tiempo 5-23, 25- 
33 y 33-53 s, en el intervalo de 25 a 33 s se observa cómo se libera la mayor parte 
de energía que produce la ruptura de la segunda aspereza; d) mapa de las 
diferentes estaciones utilizadas en el estudio; e) diagramas tiempo-deslizamiento 
en las diferentes estaciones, en negro los observados y rojo los modelados. 
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Figura 2 Evolución de la ruptura en un espacio de tiempo y liberación de la energía sísmica 

Fuente: Noquet et al. (2016) 


En este contexto es necesario realizar prospecciones geotécnicas, debido 
a que en su mayoría las estructuras fallaron por efecto del suelo, por no poder 
definir un adecuado factor de amplificación del suelo, de ahí la necesidad de 
realizar estudios con interacción suelo estructura en las zonas anteriormente 
descritas. 


Este estudio trabaja con un edificio de doce pisos ubicado en la ciudad de 
Manta, al cual se realiza un análisis sísmico en coordenadas de piso, se 
consideran tres grados de libertad por planta y se emplea el espectro de Manta en 
componentes E-W y N-S, este primer modelo es con base empotrada y determina 
los desplazamientos y fuerzas en el centro de masas. 

Se presenta el programa desp_fuer_NP que determina los 
desplazamientos y fuerzas de cada pórtico que conforma la estructura. Para 
obtener el desplazamiento y fuerza resultante se aplicó el Criterio de Combinación 
Modal de la Norma Técnica de Perú de 2003. (Agular y Ortiz, 2016). 


254 


Roberto Aguiar, Dennis Ortiz 


ASCE-10 propone mediante ecuaciones obtener un modelo de un grado 
de libertad equivalente de cualquier estructura; de esta modo se puede estimar la 
rigidez, masa, altura, período, desplazamiento y fuerza equivalente. 
Posteriormente al sistema de Igdl equivalente se lo hace interactuar con el suelo; 
a través de la guía NEHRP-12 se detalla de manera minuciosa el cálculo de los 
factores de rigidez y amortiguamiento que se toman en cuenta para el 
modelamiento de una estructura con vigas de cimentación. 

Se presenta también el programa rig¡dez_suelo , que determina la rigidez 
del suelo y la velocidad de onda de corte; y de igual forma el programa 
resortes_ZS, que calcula los resortes de la cimentación para realizar interacción 
suelo estructura. (Aguiar y Gallegos, 2014). 

Por último se pretende relacionar el nivel de daño que sufrió la edificación 
que se indica en la figura 3, a partir de las derivas de piso calculadas en la 
estructura. Se encerró con círculo rojo las diferentes afectaciones registradas. 



e) f) 

Figura 3 a) Vista posterior; b) Vista Lateral; c), d) Fachada principal con entrantes; e), f) Daños en 

manipostería 
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2. MÉTODO ESPECTRAL EVALUADO EN CENTRO DE MASAS 


El Análisis Modal Espectral es uno de los métodos más utilizados para el 
análisis sísmico de estructuras, en la práctica profesional para el rango elástico, el 
método en mención es muy eficaz y con resultados precisos. Se detalla el método 
en forma general, aplicable a cualquier estructura. (Aguiar, 2012). 

En primer lugar, es imprescindible calcular la matriz de rigidez lateral de 
todos los pórticos que conforman la estructura, a continuación se calcula la matriz 
en coordenadas de piso, para proceder con el cálculo de la matriz de masas la 
cual está en función de las cargas de la estructura, con estos dos últimos 
parámetros deducidos se obtienen las propiedades dinámicas y modos de 
vibración de la estructura los mismos que dan paso para determinar los factores 
de participación modal y las aceleraciones espectrales, en combinación de estos 
elementos calculados se obtienen los desplazamientos máximos modales y las 
fuerzas máximas modales. 


2.1. Matriz de Rigidez en Coordenadas de Piso 

La matriz de rigidez se calcula con un modelo de losa rígida con tres 
grados de libertad por planta considerados en el centro de masas (CM), dos 
componentes de desplazamiento según los ejes X, Y, y una rotación con respecto 
a un eje perpendicular a la losa. 


La matriz de rigidez en coordenadas de piso, es simétrica, razón por la 
cual solo se precisa calcular las sub matrices de la matriz triangular superior. 


k e 


K xx 


Kxy 

Kyy 


Kxe 

Kye 

Kee 


( 1 ) 


2.2. Matriz de Masas en Coordenadas de Piso 


La matriz de masas de cualquier estructura se halla a partir de la 
evaluación de la energía cinética, T c en el centro de masas. 

Para el modelo a analizar, la matriz de masas es la siguiente. 


~m 1 


\h 

m 2 


h 


j = 


m t 


Ji 

m n . 


Jn- 


Donde la masa total del piso i; m n la masa del último piso; J t es el 
momento de inercia de la masa m ¿ ; ) n es el momento de inercia de la masa m n . Se 
destaca que son momentos de inercia con respecto al centro de masa de cada 
piso. Para el piso ‘i’ se tiene las ecuaciones (3) y (4). 
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m.i = 


C Tj * A t 

9 


:3) 


Donde CT t es la carga total que actúa en ese piso; A í es el área del piso, y 


m¡ , , 

J i ~ 12 ( a ‘ + ) 


; 4 ) 


Donde m¡ es la masa del piso ¿; a¿, son las dimensiones de la losa en el 


piso i. 


2.3. Desplazamientos Máximos 

Para obtener los desplazamientos inelásticos en el centro de masas de 
una estructura en cada modo de vibración, se utiliza la ec. (5). 

q(í) = r(O © 2j4di0(O (5) 

Donde r® es el factor de participación modal en valor absoluto y 
representa que tanto participa el modo en la respuesta; es el período de 
vibración del modo ¿; A di es la aceleración espectral asociada al período T ¡; 0® es 
el modo de vibración en el modo i. Para hallar la respuesta total se debe utilizar un 
criterio de combinación modal. 

2.4. Fuerzas Máximas Modales 

Para encontrar las fuerzas en coordenadas de piso en cada modo de 
vibración se emplea la ec. (6). 

Q(i) = r^A di M 0® ( 6 ) 

Al igual que con los desplazamientos para hallar la respuesta total se debe 
aplicar un criterio de combinación modal. Cabe destacar que para hallar las 
fuerzas, las aceleraciones previamente deben ser afectadas por el factor de 
reducción de las fuerzas sísmicas R. 

2.5. Criterios De Combinación Modal CQC 

El criterio CQC (Complete Quadratic Combination), Chopra (2001), 
considera la posibilidad de acoplamiento entre los modos de vibración. 

N N 

r2= ZZ Pt ' r ‘ r ' (7) 

i= 1 7=1 

2 (1 + a)a 15 _W nj 

(1 - a 2 ) 2 + 4f 2 a(l + a) 2 a ~ W ni 


Pij = 
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= BjhTjfa + 

Pii (1 - a 2 ) 2 + 4fája(l + a 2 ) + 4(f ¿ 2 + f/)a 2 

Se presentan dos ecuaciones para p ¡; , la primera considera un factor de 
amortiguamiento igual en todos los modos de vibración. Y la segunda considera 
factores de amortiguamiento diferentes en cada modo de vibración. La validez de 
las dos ecuaciones dependerá de los datos de factor de amortiguamiento con los 
que cuente el investigador. 

Cuando las frecuencias están bastante separadas, el criterio de la 
combinación cuadrática completa, proporciona valores similares al criterio del 
máximo valor probable. 

Para aplicar el criterio de combinación modal cuando se realice un análisis 
en coordenadas de piso, se debe dividir el número de modos de vibración para 
tres y de esa manera poder calcular por separado las fuerzas para el sentido X, Y 
y rotación. 

3. MÉTODO ESPECTRAL A NIVEL DE PÓRTICO 

A partir de los desplazamientos en cada modo de vibración se realiza un 
análisis a nivel de pórtico, para los desplazamientos en cada pórtico se utiliza la 
ec. (8). 


pd) = ¿4® * q(0 ( 8 ) 

Donde A t es la matriz de compatibilidad del pórtico i, que relaciona las 
coordenadas laterales de un pórtico con las coordenadas de piso de la estructura. 
(Agular 2008). La forma de la matriz A (i) es la siguiente: 

'cosa sena 

^ (í) _ cosa sena 

cosa sena r n J 

Donde a es el ángulo que forma la orientación positiva del pórtico con el 
eje de las X. Para pórticos paralelos al eje X este ángulo vale 0 o y para pórticos 
perpendiculares al eje X vale 90°. Por otra parte r n es la distancia desde el centro 
de masas al pórtico en el piso n, será positiva si la orientación del pórtico rota con 
respecto al centro de masa en sentido antihorario. La matriz A tiene np filas y 3 * 
np columnas, donde np es el número de pisos del pórtico. 

Los desplazamientos que se obtienen son en cada modo de vibración; se 
aplica el criterio de combinación modal para obtener el desplazamiento resultante 
en cada altura de piso. Para las fuerzas a nivel de pórtico se utiliza la ec. (10). 


r i 


r 2 


| 9 ) 


Pi = KL i * Vi 


( 10 ) 



258 


Roberto Aguiar, Dennis Ortiz 


Donde KL t es la matriz de rigidez lateral del pórtico i; para la fuerza 
resultante para cada altura es indispensable aplicar un criterio de combinación 
modal. 

3.1. Norma Técnica De Perú 2003 

En la Norma Técnica de Perú de 2003 se combinan los resultados 
obtenidos, en cada uno de los modos de vibración, con la siguiente ecuación: 



( 11 ) 


En la Normativa de Perú se reconoce que el criterio del valor máximo 
probable reporta valores bajos y que el criterio de superposición directa da valores 
muy altos por lo que lo más conveniente es combinar estos dos criterios en forma 
lineal con los coeficientes indicados en la ecuación. 

Programa CEINCI-LAB 

Se describe el programa de CEINCI-LAB, que se utilizó para obtener los 
desplazamientos y fuerzas de cada pórtico de la estructura. (Aguiar y Ortiz, 2016). 

El programa desp_fuer_NP, devuelve un vector compuesto por. 
p Desplazamientos en el pórtico por cada modo de vibración. 


pt Desplazamientos resultantes luego de la combinación modal. 

P Fuerzas en el pórtico por cada modo de vibración. 

PT Fuerzas resultantes luego de la combinación modal. 


A partir de los siguientes datos de entrada. 


KL Matriz que contiene las matrices de rigidez de cada pórtico. 

A Matriz que contiene las matrices A de cada pórtico. 

np Número de pisos de la estructura. 

q Es la matriz que contiene los desplazamientos en cada modo. 

Npor Indica el pórtico en el que realiza el cálculo. 


4. MODELO SIMPLIFICADO DEL ASCE-10 

El Modelo Simplificado es un sistema de un grado de libertad el mismo 
que es equivalente al de la estructura de múltiples grados de libertad. Es en este 
sistema que se encuentra la respuesta dinámica ante una acción sísmica dada. 

Existen varios modelos para definir el sistema equivalente de Igdl, todos 
ellos se derivan de los principios básicos de la dinámica de estructuras los 
modelos que aquí se indican son los propuestos por; Rodríguez (1994) del 
Instituto de Ingeniería de la UNAM, México; Fajfar y Gaspersic (1996) de la 
Universidad de Ljubljana, Slovenia; Esteva (1999) del Instituto de Ingeniería de la 
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UNAM, México; Aguiar (2001) del Centro de Investigaciones Científicas de la 
ESPE, Ecuador y Ayala (2001) del Instituto de Ingeniería de la UNAM, México. 

El modelo utilizado por la ASCE-10, utiliza los períodos y el primer modo 
de vibración en función de Kxx que corresponde a la sumatoria de las matrices de 
rigidez lateral en el sentido de análisis y m matriz de masas de la estructura. Para 
transformar una estructura de múltiples grados de libertad a una estructura de un 
grado de libertad se emplean las siguientes ecuaciones: 


hj 


m12. 

mili 

mIOí 

m9 

mG 

iíi7 

mG 

m5 

m4 

m3 

rn2 

mi 


\2^ 
|1 . 


meqi 


he 


Ceq 

Keq 


Figura 4 Modelo de múltiples grados de libertad a un grado de libertad equivalente 
Masa Equivalente 


M ef = 


[0' efl *M *7'] 2 

tyefl * M * 0e/l] 


Rigidez Equivalente 


Kef - M e f * 



Altura Equivalente 


_ [0' efl * M * h'] 
ef ~ [0' efl *M *J'] 


( 12 ) 


(13) 


(14) 


Frecuencia Equivalente 


Wn ef = 



(15) 
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Período Equivalente 

2n 

To = w^r f 

Desplazamiento Equivalente 

(T 0 \ 2 

q ° = \2ÜJ * Mef 

Fuerza Equivalente 

F 0 F e f * c [ 0 


(16) 


5. CÁLCULO DE RIGIDEZ Y AMORTIGUAMIENTO DEL SUELO 
5.1. Rigidez Del Suelo 

En diversas recomendaciones Internacionales, tales como las del NEHRP, 
la calidad del suelo superficial se indica mediante el promedio de la velocidad de 
propagación de las ondas de corte en los primeros 30 m de profundidad Fs30. De 
manera similar, la norma ecuatoriana de la construcción, en el capítulo de peligro 
sísmico NEC-15, proporciona la tabla 1 donde utiliza la velocidad promedio de las 
ondas de corte para la definición de los perfiles de suelo. 


Tabla 1 Criterios para clasificar suelos 
Fuente: NEC-15, peligro sísmico 


Tipo de Perfil 

v s 

NoN rh 

s„ 

C 

entre 360 y 760 m/s 

mayor que 50 

mayor que 100 kPa 

D 

entre 180 y 360 m/s 

entre 15 y 50 

entre 100 y 50 kPa 

E 

menor 180 m/s 

menor de 15 

menor de 50 kPa 


Tabla 2 Clasificación sísmica de suelos 


Fuente: Código Internacional de Construcción 


Tipo 

Vs3o (m/s) 

Descripción 

A 

>1,500 

Roca dura 

B 

760- 1,500 

Roca 

C 

360 - 760 

Roca blanda y suelo muy denso 

D 

180-360 

Suelo rígido 

E 

<180 

Suelo blando 

F 

Suelos especiales (licuables, colapsables, arcillas 
muy plásticas, suelos orgánicos) 


El módulo de rigidez a corte del suelo es el parámetro determinante en las 
propiedades del suelo, un elemento diferencial del suelo pierde rigidez conforme la 
onda sísmica sube a la superficie, es decir conforme la distorsión angular del 
elemento aumenta, la rigidez a corte del suelo disminuye. Ver figura 5a. 
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Por lo cual el amortiguamiento aumenta, conforme el módulo de corte 
disminuye, problema del suelo de México. Ver figura 5b. (Gazetas, 201 2). 



Figura 5 a) Módulo de rigidez a corte Gmax; b) Degradación de la rigidez del suelo 
Fuente: Gazetas (2012) 


Para obtener estos parámetros se debería realizar ensayos de laboratorio, 
la ASCE 10 simplifico esto en función del inicio del espectro conocido como 
S DS / 2.5 y presenta la tabla 19.2-1, para efectos de este estudio se realizó una 
adaptación de la tabla mencionada y se indica la tabla 3 la cual a partir de la 
aceleración máxima en suelo PGA permite contemplar la degradación del módulo 
de rigidez G. Para valores medios se debe interpolar y para suelo tipo F es 
necesario realizar ensayos. 


Tabla 3 Factor de reducción de rigidez del suelo 


Tipo de 
Suelo 

G/Gmax 

PGA 


<0.1 

0.4 

>0.8 

A 

1,00 

1,00 

1,00 

B 

1,00 

0,95 

0,90 

C 

0,95 

0,75 

0,60 

D 

0,90 

0,50 

0,10 

E 

0,60 

0,05 

a 

F 

a 

a 

a 


El módulo de rigidez a corte inicial o máximo G max está en función de la 
velocidad de onda de corte V s 30, parámetro dado como dato a partir de definir el 
tipo de suelo. 


Gmax = r(Vs30y/g (19) 

Donde Vs30 velocidad de onda de corte promedio a 30 m de 
profundidad; / peso específico del suelo; g gravedad. 
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Se obtiene la relación G/G max , y con esa proporción se calcula el módulo 
de corte del suelo con degradación de sismo G. 

G = % G max (20) 

Por último la velocidad de onda de corte con degradación de sismo Vs. 

Vs = V(G *g)/r ( 21 ) 

5.2. Resortes Para Cimentación 

La guía NEHRP-NIST 2012 sugiere aplicar las fórmulas de Gazetas (1991) 
y modificadas por Mylonakis (2006), estas fórmulas se orientan de modo que L > 
B, es decir que la rigidez en x siempre será en dirección de la dirección mayor. Ver 
la figura 6. 



Figura 6 Configuración de parámetros para resortes 
Fuente: Gazetas (2012) 


La rigidez del suelo con interacción del suelo suele definirse como 
funciones de impedancia que representa la frecuencia dependiente de la rigidez y 
el amortiguamiento característico de la interacción suelo estructura. 


Kj = kj + iwXj ( 22 ) 

Pero la NERPH (NIST2012) separa la rigidez del amortiguamiento, así la 
rigidez depende de las características de la cimentación: 


k = K r] a ( 23 ) 

Donde K rigidez superficial de la cimentación; r¡ factor que toma en cuenta 
el desplante; a factor de rigidez dinámica que toma en cuenta la vibración de la 
estructura. 

La tabla 4 calcula los resortes en función de las características de la 
cimentación. Se destaca que existe una tabla similar para los amortiguamientos. 
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Tabla 4 Rigideces 

para zapatas aisladas y vigas de cimentación. 

SENTIDO 

RIGIDEZ SUPERFICIAL 

FACTOR DE DESPLANTE 

FACTOR RIGIDEZ DINÁMICA 

Rigidez z 

K 2GL Í0 73 + 1 54Í' B 'l°' 75 l 

Vz = [l 

+ JL( 1 + 13 B )lf 1 + 0 2(^f 


a z = 1- 

(0.4 + ^|)a o 2 


K * ~l-v [° 73 + I' 54 UJ J 

+ 21Bl 1 + lá ijJ 1 + U2 U Bl) 

[(l + 3(L/5- l)) + a ° 

Rigidez y 

2 GL [ /fí\°- 85 l 

^ = _l 2+2 . 5 (_) j 

Vy = 

LD 

1 + 0.15 ñ 
B 



a Y = 1 

Rigidez x 

K _ K °- 2GL h ( B )l 

Vx = 

D 

1 + 0 - 15 J B 

f 1 + 0 52(^ri 


a x = l 

Kx ~ Ky 0.75- vL 1 UJJ 

[ 1 + Ub¿ \BL 7 ) \ 

Torsión z 

r / B\ 10 i 

K ZZ =GJ 0 - 7S 4+11(1--) 

— K)(^f 

a zz ~ 1 

(0.33 - 0.03 JL/B - 1 )cl 0 2 
( 0.8 \ ? 
\l + 0.33 (L/B - 1)J + U ° \ 


Rotación 

yy 

*"=1 ^°- 75 [ 3 (¿n 

/d w \°- 6 r /d w \ 19 

= 1 + 0-92 (f) |l.5 + (f) 

1 

a YY = l- 

0.55 a 0 2 

[(°' 6+ (¿/B) 3 ) + a » 2 J 


Rotación 

XX 


r, xx = 1 + 1.26^[ 



«xx = 1 - 

(0. 55 + O.Ol^/L/5 — l)a 0 2 

( 2 - 4+ a/By ) +a ° 2 \ 



Donde G módulo de rigidez degradada del suelo; a Q factor adimensional 
que toma en cuenta la frecuencia de la estructura. 


Donde w frecuencia de vibración de la estructura generalmente la del 
primer modo de vibración; Vs velocidad de onda de corte con efecto de 
degradación de la rigidez del suelo. Las restantes variables de la tabla 4 están 
definidas en forma gráfica en la figura 6. 

Programas CEINCI-LAB 

Se describen los programas de CEINCI-LAB, que se utilizaron para 
obtener la rigidez de los resortes que se aplican en la interacción suelo estructura. 
(Aguiar y Gallegos, 2014). 

En primer lugar se utilizó el programa rigidez_suelo, el cual devuelve un 
vector compuesto por. 

G Módulo de corte del suelo con degradación de sismo en t/m 2 . 

v Peso específico del suelo en t/m 3 . 

Vs Velocidad de onda de corte con degradación de sismo m/s. 

A partir del ingreso los siguientes datos de entrada. 
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„ Vector que contiene tipo de suelo. 

A=1; B=2; C=3; D=4; E=5. 

V s 30 Velocidad de onda de corte Vs30 en m/s 
p Peso específico del suelo en t/m 3 . 

PGA Aceleración máxima como una fracción de la gravedad 
V Módulo de poisson del suelo. 

Posteriormente se utiliza la función resortes_ZS, la misma que otorga un 
vector compuesto por. 

Kd Rigidez horizontal en t/m . 

Kr Rigidez rotacional en t * m . 

Cd Amortiguamiento horizontal en t * s/m. 

Cr Amortiguamiento rotacional 

Con los siguientes datos de entrada. 

Z Vector que contiene características de la viga de cimentación. 

Viga de Cimentación 

LT Longitud en el sentido de análisis en m eje X. 

AT Long. transversal (ancho de zapata) al sentido de análisis en m. 

h Altura total de la viga en m. 
e Espesor de la zapata en m. 
b Ancho de la viga en m. 

D Altura de desplante de la viga de cimentación en m. 

G Vector que contiene G, v, Vs. 

OM Frecuencias modales. 


6 . MODELO DE INTERACCIÓN SUELO ESTRUCTURA 

Para condiciones de suelo empotrados la rigidez del suelo es infinito, de 
este modo se tiene la ecuación (25). 


q = F/K 


(25) 



Figura 7 Modelo de cuerpo libre con efecto de cimentación. 
Fuente: Gallegos (2014) 
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Teniendo en cuenta que para sistemas de un grado de libertad se tiene, 
(Clough and Penzein, 1993) ver la figura 7: 


_ (i 

Jm 

(26) 

2n ira 

T = — = 2n — 
a) \ k 

(27) 

T 2 = ( 2n) 2 - = (2n) 2 -¿- 
F/q F 

(28) 

Considerando interacción suelo estructura el 
determina por la ecuación (29). 

desplazamiento final se 

Qs Qo T q d + q r 

(29) 

Donde q 0 desplazamiento de la estructura con base empotrada; q d es el 
desplazamiento horizontal del suelo, q r es el desplazamiento rotacional del suelo. 

II 

^1? 

(30) 

F a h 

c¡ r — ah — —— — h 

K r 

(31) 


h = h e f + D 


Donde K d rigidez de desplazamiento horizontal de la cimentación con 
interacción del suelo; K r rigidez de rotación de la cimentación con interacción del 
suelo; h es la sumatoria de la altura equivalente del sistema de un grado de 
libertad h ef más la altura de desplante D. 

Con interacción suelo estructura el período final, de manera aproximada 
está definido por la ecuación (32). (Aguiar, 2012). 

T 2 = T 0 2 + T d 2 + T 2 ( 32 ) 

Donde T 0 período de vibración la estructura con base empotrada; T d es el 
período de vibración traslacional considerando a la estructura infinitamente rígida 
con un solo grado de libertad; T r es el período de vibración rotacional 
considerando a la estructura totalmente rígida con un solo grado de libertad. 



T d = 2n * 


(33) 
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(34) 


Donde m t es la masa total del sistema; J es el momento de inercia 
rotacional con respecto a la base de la cimentación; J = T l m i h i 2 . Las restantes 
variables han sido ya definidas. 

La ecuación (32) se obtiene a partir de la fórmula propuesta por la ASCE 
7-10 en el capítulo 19 de interacción suelo estructura para diseño sísmico, la 
misma se encuentra definida por: 


T s I k kh 2 

V 0 = j 1 + r d + l^ 


(35) 


7 . ANÁLISIS DE LA ESTRUCTURA 

Se analiza un edificio de 12 pisos de hormigón armado, cimentada sobre 
vigas de cimentación, módulo de elasticidad del hormigón E = 15000 * /f'c y 
fe = 210 kg/cm 2 , ubicado en la ciudad de Manta. 

En la figura 8 se indican las plantas que representan al edificio. Se 
considera sentido transversal a aquel que está en dirección de los pórticos 
identificados por letras y, sentido longitudinal, al ortogonal es decir al que está 
identificado por números. La estructura es simétrica en el sentido longitudinal y se 
conforma por muros de corte en el ducto de gradas, ascensor y parte del eje H. 
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Figura 8 a) Plantas Pisos 1 , 2 y 3 - Tipo I; b) Planta Piso 4 - Tipo II; c) Plantas Pisos 5, 6, 7, 8, 9 y 
1 0 - Tipo 1 1 1 ; d) Planta Piso 1 1 - Tipo IV; e) Planta Piso 1 2 - Tipo V 


La carga por metro cuadrado de la losa de 20 cm, alivianada en dos 
sentidos es de 288.00 kg/m 2 , la carga de manipostería, vigas, columnas y muros 
varía según el tipo de losa, la carga que aportan los acabados es de 60.00 kg/m 2 ; 
de esta forma se obtuvieron los valores para carga muerta, la carga viva es de 
200.00 kg/m 2 . 

Para el análisis sísmico se considera el 25% de la carga viva, por la 
probabilidad de ocurrencia de un terremoto con toda la carga viva. Se determina la 
masa con la ec. (3) y momento de inercia de la masa con la ec. (4). Así la tabla 5 
refleja los datos con los cuales se efectuó el análisis. 


Tabla 5 Masas e Inercias de la estructura. 


Losa 

D ( t/m 2 ) 

L (t/m 2 ) 

T (t/m 2 ) 

Área (m 2 ) 

m (ts 2 /ni) 

J ( tms ) 

Tipo 1 

1,060 

0,200 

1,110 

446,780 

50,605 

3768,193 

Tipo II 

1,000 

0,200 

1,050 

731,640 

78,390 

9558,877 

Tipo III 

0,960 

0,200 

1,010 

574,970 

59,257 

5678,510 

Tipo IV 

0,850 

0,200 

0,900 

376,460 

34,573 

2169,216 

TipoV 

0,400 

0,200 

0,450 

18,430 

0,846 

2,599 
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En la figura 9 se presenta de manera indistinta pórticos que forman parte 
de la estructura. Para el sentido transversal los pórticos E y D, y para el sentido 
longitudinal los pórticos 5, 2 y 9. 


(j) (*) <p C&) Cp <p Cp <íp 



Figura 9 a) Pórtico E; b) Pórtico D; c) Pórtico 5; d) Pórtico 2 y 9 


7.1. Análisis Estructura Con Base Empotrada 

Con el propósito de que la matriz de masas sea diagonal, las coordenadas 
de piso se consideraron en el centro de masas. Se calculó la matriz de rigidez con 
un modelo de losa rígida con tres grados de libertad por planta, dos componentes 
de desplazamiento según los ejes X, Y, y una rotación con respecto a un eje 
perpendicular a la losa. Las sub matrices de la matriz de masas son: 



M = 
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50.60 


50.60 


50.60 


78.39 


59.26 


59.26 


59.26 


59.26 


59.26 


59.26 


34.57 


0.85 


3768.19 


9558.88 


5678.51 


5678.51 


5678.51 


5678.51 


5678.51 


5678.51 


2.60 


Propiedades Dinámicas 

En la tabla 6, se indican los períodos de vibración de la estructura con 
base empotrada. El período fundamental es 1.155 s. 


Tabla 6 Períodos de vibración de la estructura con base empotrada 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1,155 

1,138 

0,888 

0,381 

0,336 

0,286 

0,201 

0,181 

0,164 

0,162 

0,161 

0,119 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

0,107 

0,101 

0,093 

0,085 

0,078 

0,073 

0,064 

0,063 

0,058 

0,053 

0,051 

0,050 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

0,045 

0,043 

0,042 

0,040 

0,039 

0,037 

0,037 

0,035 

0,034 

0,033 

0,031 

0,030 


Espectro de Análisis 

Se muestra los espectros que se obtuvieron en la ciudad de Manta para el 
terremoto del 16 de abril de 2016, se aprecia la presencia de picos característicos, 
típicos de los sismos interplaca tipo thrust. Ver figura 1 0. (Aguiar 201 0). 


Espectros Manta 16-04-2016 



Figura 10 Espectros para la dudad de Manta, terremoto del 16 de abril de 2016. 
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Los picos característicos son aquellos valores en los cuales la aceleración 
espectral vuelve a subir y por ende ingresa una mayor cantidad de energía a la 
estructura. 


Método Espectral 

Se realiza el análisis sísmico para las dos componentes del terremoto del 
16 de abril de 2016, la condición más desfavorable corresponde a la componente 
E-W. 


Se trabaja con un factor de reducción de las fuerzas sísmicas R = 4 esto 
se debe a que se considera un factor de reducción de resistencia por ductilidad 
= 2, factor de sobre resistencia R^ = 2, y factor de redundancia R r = 1. 

Las respuestas máxima probables en el centro de masas de cada uno de 
los piso se indican en la tabla 7. Donde q es el desplazamiento inelástico, F es la 
fuerza lateral y y es la deriva de piso. 


Tabla 7 Respuestas Máximas Probables 



Análisis Transversal 

Análisis Longitudinal 


COMPONENTE E-W 

COMPONENTE N-S 

COMPONENTE E-W 

COMPONENTE N-S 

piso 

q (cm) 

F(T) 

Y (%) 

q (cm) 

F(T) 

Y (%) 

q (cm) 

F(T) 

Y (%) 

q (cm) 

F(T) 

Y (%) 

1 

0,93 

47,21 

0,29 

0,85 

49,09 

0,27 

1,91 

57,99 

0,60 

1,54 

49,62 

0,48 

2 

2,56 

71,34 

0,51 

2,29 

79,59 

0,45 

3,62 

94,38 

0,54 

2,92 

77,89 

0,43 

3 

4,35 

71,40 

0,56 

3,76 

90,43 

0,46 

3,78 

79,13 

0,05 

3,05 

67,60 

0,04 

4 

4,98 

113,09 

0,20 

4,11 

127,17 

0,11 

5,52 

131,80 

0,54 

4,43 

115,14 

0,43 

5 

6,98 

89,85 

0,62 

5,58 

97,96 

0,46 

7,68 

101,10 

0,67 

6,12 

90,45 

0,53 

6 

9,03 

93,67 

0,64 

7,07 

97,91 

0,47 

9,94 

107,90 

0,71 

7,88 

92,22 

0,55 

7 

10,99 

95,13 

0,61 

8,47 

91,37 

0,44 

12,15 

104,87 

0,69 

9,59 

85,12 

0,53 

8 

12,74 

96,39 

0,55 

9,71 

83,84 

0,39 

14,24 

97,04 

0,65 

11,21 

76,58 

0,51 

9 

14,43 

98,26 

0,53 

11,05 

86,39 

0,42 

16,41 

95,95 

0,68 

12,93 

78,30 

0,54 

10 

16,09 

108,30 

0,52 

12,45 

103,31 

0,44 

18,68 

130,21 

0,71 

14,75 

108,58 

0,57 

11 

19,57 

97,12 

1,09 

15,06 

102,38 

0,81 

20,92 

110,53 

0,70 

16,54 

92,07 

0,56 

12 

24,05 

11,66 

1,40 

18,82 

9,02 

1,18 

23,87 

8,90 

0,92 

19,00 

7,47 

0,77 


Mayores derivas de piso se dan en la planta 12 debido a que existen solo 
cuatro columnas; este piso no se contempló en el diseño original. Ver figura 8e. 

Con derivas de piso inelásticas mayores a 0.1% se presenta daño en la 
mampostería (Filiatrault, 2015). Las derivas de piso encontradas justifican el daño 
que se presentó en la mampostería de la estructura. 

En la figura 1 1 se presentan las derivas de piso encontradas en el análisis 
sísmico en sentido transversal y longitudinal. 
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Figura 11 Derivas de piso: a) Sentido Transversal E-W; b) Sentido Transversal N-S; 
c) Sentido Longitudinal E-W; d) Sentido Longitudinal N-S. 


7.2. Método Espectral a Nivel de Pórtico 

En la tabla 8 se presenta la respuesta dinámica espectral tanto del pórtico 
D, ver figura 9b, como del pórtico 2, ver figura 9d. Donde p es el desplazamiento y 
P es la fuerza lateral. 

El pórtico D es de siete pisos razón por la cual en los niveles 1 , 2, 3, 1 1 y 
12 el valor es cero, esto ocurre también en el pórtico 2 en los niveles 11 y 12. 

En comparación con la tabla 7 donde se obtienen los desplazamientos y 
fuerzas en el centro de masas, se puede apreciar que los desplazamientos a nivel 
de pórtico son similares y las fuerzas son proporcionales al número de pórticos 
que existe en el sentido de análisis. 
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Tabla 8 Desplazamientos y Fuerzas en el Pórtico D y Pórtico 2 



Análisis Sentido Transversal 

Análisis Sentido Longitudinal 

Pórtico D 

Pórtico 2 

COMPONENTE E-W 

COMPONENTE N-S 

COMPONENTE E-W 

COMPONENTE N-S 

piso 

P (cm) 

P(T) 

P ( cm ) 

P(T) 

P (cm) 

P(T) 

P (cm) 

P (T) 

1 

0,00 

0,00 

0,00 

0,00 

1,84 

5,27 

1,50 

4,51 

2 

0,00 

0,00 

0,00 

0,00 

3,49 

8,58 

2,82 

7,08 

3 

0,00 

0,00 

0,00 

0,00 

3,62 

7,19 

2,92 

6,15 

4 

2,29 

18,85 

1,97 

21,20 

5,30 

11,98 

4,25 

10,47 

5 

3,82 

14,98 

3,45 

16,33 

7,43 

11,23 

5,93 

10,05 

6 

5,39 

15,61 

5,02 

16,32 

9,65 

11,99 

7,66 

10,25 

7 

6,96 

15,85 

6,59 

15,23 

11,84 

11,65 

9,35 

9,46 

8 

8,33 

16,06 

7,98 

13,97 

13,91 

10,78 

10,96 

8,51 

9 

9,62 

16,38 

9,26 

14,40 

16,07 

10,66 

12,68 

8,70 

10 

11,09 

18,05 

10,63 

17,22 

18,34 

14,47 

14,50 

12,06 

11 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

12 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 

0,00 


7.3. Sistema de Un Grado de Libertad Equivalente 

Para determinar los parámetros de un sistema de un grado de libertad 
equivalente se precisa obtener los períodos y el primer modo de vibración 
equivalente, mismos que se representan en la tabla 9, para lo cual se trabaja con 
la sub matriz Kxx de la ec. (1) y la matriz de masas m de los dos sentidos de 
análisis, motivo por el cual se obtiene valores diferentes. 


Tabla 9 Períodos y Modos de Vibración Equivalentes 



Análisis Sentido Transversal 

Análisis Sentido Long 

itudinal 

Períodos 

i 

2 

3 

4 

5 

6 

i 

2 

3 

4 

5 

6 

1,061 

0,342 

0,189 

0,160 

0,118 

0,085 

1,138 

0,336 

0,181 

0,163 

0,101 

0,073 

7 

8 

9 

10 

11 

12 

7 

8 

9 

10 

11 

12 

0,065 

0,052 

0,045 

0,039 

0,035 

0,034 

0,058 

0,050 

0,043 

0,037 

0,034 

0,030 

Primer Modo 
de Vibración 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

-0,002 

-0,005 

-0,009 

-0,015 

-0,023 

-0,033 

0,006 

0,011 

0,011 

0,017 

0,025 

0,033 

7 

8 

9 

10 

11 

12 

7 

8 

9 

10 

11 

12 

-0,042 

-0,050 

-0,058 

-0,065 

-0,071 

-0,075 

0,041 

0,049 

0,057 

0,064 

0,071 

0,076 


Se determina la tabla 10 a partir de emplear la ecs. (12-18), se puede 
apreciar que la altura equivalente h ef de los dos sentidos de análisis en promedio 
es de 25,00 m. que corresponde de manera proporcional al octavo nivel de la 
estructura; al comparar con la tabla 7 se observa que los desplazamientos en el 
piso 8 son muy semejantes a los que se calcularon a través del sistema de un 
grado de libertad equivalente. 
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Tabla 10 Sistema de Un Grado de Li 

bertad Equivalente 


A. S. Transversal 

A. S. Longitudinal 


COM E-W 

COM N-S 

COM E-W 

COM N-S 

K ef (t/m) 

14759,05 

13796,47 

M ef (ts 2 /m) 

420,63 

452,97 

h ef (m) 

25,50 

24,72 

Wn ef (Hz) 

5,92 

5,52 

T„(s) 

1,06 

1,14 

q« (cm) 

12,71 

10,74 

13,72 

10,80 

F„(T) 

1876,00 

1584,97 

1892,83 

1490,65 


8. Análisis Estructural Con Interacción Suelo Estructura 


Para la ciudad de Manta el tipo de suelo que corresponde es C, la 
velocidad de onda de corte Vs30 será igual a la que proporciona la NEC-15, en el 
capítulo de peligro sísmico, peso específico del suelo y = 1.60 t/m 3 , coeficiente de 
poisson v = 0.28, la aceleración inicial del espectro elástico PGA es igual al 
máximo valor del registro sísmico de Manta fraccionado para la gravedad g 
expresada en gals, en las componentes E-W y N-S. 

La tabla 11 determina el módulo de rigidez y la velocidad de onda de corte 
del suelo, para lo cual se aplicó las ecs. (19-21). 


Tabla 11 Rigidez del suelo para un perfil tipo C, en componentes E-W y N-S 


Tipo de Perfil 

Vs30 

Componente 

PGA 

Gmax 

G/G max 

G 

Vs 

C 

360 

E-W 

0,4052 

21159,18 

0,748 

15827,94 

311,36 

N-S 

0,5248 

0,703 

14879,47 

301,89 


Rigideces de la Cimentación de la Estructura 

En la figura 12 se presentan las características y configuración de la 
cimentación. 



a) 
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Figura 12 a) Corte Longitudinal De Cimentación En Eje I; b) Planta de Cimentación N-14.40 


Tanto la rigidez horizontal como la rigidez rotacional dependen de los 
parámetros que proporciona la tabla 2, para ilustrar cómo se realiza este estudio 
se presenta de manera particular el cálculo de la viga de cimentación del eje I, en 
componente E-W. 


Para rotación alrededor del eje Y de la viga de cimentación figura 12. 


LT = Longitud total de la cimentación = 18.10 m, L = LT/2 = 9.05 m, AT= ancho 
total de la cimentación = 1 .40 m, B = A T/2 = 0.70 m, h = altura total de la viga = 1 .20 m, 
e = espesor zapata = 0.40 m, b = ancho viga = 0.60 m, D = altura de desplante viga 
cimentación = 1.60 m 


Rigidez Horizontal 

Rigidez Superficial 


2 GL 

/Bn°- 85 1 

k y = ^ 

2 + 2.5 - 

2 - v 

V L ) 


Ky = 


2 * 15827.937*9.05 


2-0.25 


2 + 2 


/0.70\°' 85 

■ 5 fes) 


373883.705 
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K y = K v ~- 


0.2GL 


K x = 373883.705 -• 


0.75 -v 
0.2 *15827.937*9.05 


M?)] 


0.75-0.25 


/0.70\l 

1 -(^05)J = 321018 ' 396 


Factor por Desplante 


Vx 


1 + 0.15 


1 + 0 


r 0 (ZwA W \ 
\ BL 2 ) 


VX = 


1 + 0.15 


1.60 


0.70 


/1.1231 * 45.52\ 0,4 
1 + 0 ' 52 ( 0.70 * 9.05 2 ) 


= 1.8361 


Factor de Rigidez Dinámica 


a x — 1 


Kd — K x * 7] x * cc x 


Kdl = 321018,396 * 1.8361 * 1 = 589426, 298 t/m 


Rigidez Rotacional 

Rigidez Superficial 


Kyy = 


K, 


YY 


1 — V 


15827.937 
1 - 0.25 


691.803 a75 



12537292.925 


Factor por Desplante 


, ^no?í dw )°' 6 íi ^í^fñ ’ 0 ' 6 

r) YY - 1 + 0.92 (— J 1-5 + (— J 


/0.7429\ 0,6 f /0.7429\ 1 ' 9 /0.70\“ a6 
*yy = 1 + 0.92 (- — ) 1-5 + (^) (^s) 


= 3.4607 
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Factor de Rigidez Dinámica 


CL 0 


wB 

~Vs 


CL 0 


5.4399* 0.70 
311.362 


0.01223 


CCyY — 1 


0.55a o 2 


(06+ 1A 1 

+ a Q 2 

U°- 6 + (l/b) 3 J 


(Xyy — 1 — 


0.55 * 0.01223 2 


^ 0.6 + 


1.4 


(9.05/0.70) 3 


+ 0.01223 2 


= 0. 99986 


Kv — Kyy * IJyY * YY 


KrI = 12537292.925 * 3.4607 * 0.9999 = 43381879. 564 t - m/rad 


La tabla 12 indica la contribución de todas las vigas de cimentación 


Tabla 12 Contribución de todas las Vigas de Cimentación 



Cimentación Sentido Transversal 

Cimentación Sentido Longitudinal 

Tipo de Perfil 

Componente 

Kd (t/m) 

Kr (t-m/rad) 

Kd (t/m) 

Kr (t-m/rad) 

C 

E-W 

4408450,55 

671280782,87 

5471274,47 

471839540,67 

N-S 

4144281,40 

631049874,40 

5143417,34 

443535346,81 


Períodos y Desplazamientos por Rigidez del Suelo 

La tabla 13 indica el aporte del suelo en el cálculo, para los períodos se 
aplicaron las ecs. (33 y 34) y para los desplazamientos las ecs. (30 y 31). 


Tabla 1 

13 Períodos y desplazamientos por Rigidez del Suelo 


A. S. Transversal 

A. S. Longitudinal 


COM E-W 

COM N-S 

COM E-W 

COM N-S 

T d (s) 

0,075 

0,067 

T r (s) 

0,058 

0,069 

q d (cm) 

0,043 

0,036 

0,035 

0,027 

q r (cm) 

0,205 

0,173 

0,278 

0,219 
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Períodos y Desplazamientos Interacción Suelo Estructura 

De este modo, la tabla 14 determina el período y desplazamiento final con 
interacción suelo estructura ec. (29) y ec. (32). 

Se puede apreciar en comparación con la tabla 10, que con interacción 
suelo estructura los períodos y desplazamientos son mayores, pero la ampliación 
por efecto de suelo es mínima. 


Tabla 14 Períodos y Despl azamientos ISE 



A. S. Transversal 

A. S. Longitudinal 

COM E-W 

COM N-S 

COM E-W 

COM N-S 

T s (s) 

1,065 

1,143 

q s (cm) 

12,959 

10,948 

14,032 

11,051 


9. CONCLUSIONES Y COMENTARIOS 

En este artículo para obtener desplazamientos y fuerzas en el centro de 
masas de la estructura se aplicó el método espectral, el mismo que dio como 
resultado derivas de piso no mayores a 0,70 % en los primeros diez pisos, lo que 
contrasta con el poco daño que sufrió la estructura; para los dos últimos niveles la 
deriva de piso excedió este valor al punto que en el nivel doce en componte E-W 
alcanzo el 1,40 %, esto se debe a la configuración de este planta la cual no se 
contempló en el diseño original y solamente posee cuatro columnas. 

Con la metodología propuesta para hallar desplazamientos y fuerzas en 
cualquier pórtico, se obtuvieron resultados similares en desplazamiento y 
proporcionales en fuerzas. Lo que afirma la validez de esta teoría. 

Finalmente en relación al modelo que consideró interacción suelo 
estructura, se apreció que el edificio de Manta por encontrarse en suelo tipo C, el 
factor de ampliación del suelo afectó muy poco a los períodos y desplazamiento, a 
su vez se debe considerar que una cimentación con Vigas T aporta un grado de 
empotramiento fuerte a la estructura, lo cual facultaría modelar la estructura sin 
tomar en cuenta la cimentación. 

El desarrollo de esta investigación ha permitido incorporar nuevos 
programas para el cálculo de la respuesta dinámica espectral de cualquier pórtico 
al sistema de computación CEINCI-LAB. Se aspira aportar con este trabajo al 
desarrollo de la ingeniería sísmica. 
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